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Abstract 


An estimate of the wavelength accuracy of a laser wavemeter is 
performed for a system consisting of a multiple-beam Fizeau interferometer 
and a linear photosensor array readout. The analysis consists of 
determining the fringe position errors which result when various noise 
sources are included in the fringe forming and detection process. Two 
methods of estimating the fringe centers are considered: 1) maximum pixel 

current location, and 2) average pixel location for two detectors with 
nearly equal output currents. Wavelength error results for these two 
methods are compared for some typical wavemeter parameters. 

Introduction 

The use of tunable laser sources in fields such as spectroscopy and 
lidar often require that some method be used to accurately determine the 
wavelength of the laser. Depending upon the constraints of the particular 
application factors such as weight, volume and spectral coverage can dictate 
what types of wavelength measuring devices are suitable to satisfy the 
experimental requirements. For those applications where space is limited 
and where both pulsed and cw lasers are used, Fabry-Perot or Fizeau 
interferometers interfaced with linear diode array sensors are attractive. 

In this approach, the wavelength is determined by measuring the fringe 
locations of the various laser sources by the diode array. Since the 
measured fringes are, in general, distorted by errors introduced in the 
detection process, the accuracy of this method depends upon how well the 
position of the measured fringes corresponds to the true fringe positions 
generated by the interferometer. 
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In this paper, the wavelength accuracy of a Fizeau-type wavemeter 
having a linear photosensor array readout is presented. The results show 
that the wavelength errors depend upon both the method used to measure the 
fringe locations and to errors in the system including: (1) response 

non-uniformity of the sensor array; (2) current quantization errors; 

(3) electronic si gnal -to-noise ratio; and (4) spectral irradiance 
non-uniformity across the interferometer aperture. 

Background 

Basically, interferometers are classified into two groups: (1) two- 

beam interferometers where the fringes can be described according to two- 
beam interference; (2) multiple-beam interferometers where fringe sharpening 
is obtained by the interference of multiple beams. A good example of a 
wavemeter of the first type is one developed by Snyder, et al. (Ref. 1) in 
which a Fizeau wedge with uncoated plates is used. This low finesse 
interferometer attains wavelength accuracy through the combined use of 
selected electronic filtering and computer algorithms to determine the 
period of interference fringes. An example of a wavemeter of the second 
type is one employing either a Fabry-Perot (Refs. 2, 3) or Fizeau interfero- 
meter (Ref. 4) in which the plates are coated for high reflectivity. This 
form of wavemeter has much sharper fringes (high finesse) than the two-beam 
type and consequently can resolve and display complex spectral structure. 

The analysis presented in this paper is for an interferometer of the 
second type. Specifically, the equations are developed for a multiple-beam 
Fizeau using a standard laser of known wavelength (e.g., frequency stabil- 
ized He-Ne) as a reference. However, the procedures used for error deter- 
mination of the fringe positions may be applied to any measurement device in 
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which the relative position of known and unknown laser fringes are measured 
via a linear photosensor array. 

Wavelength Measurement Concept 

Figure 1 shows the basic geometry of a Fizeau interferometer. In this 
diagram, the two inner surfaces of the plates are inclined at a wedge angle, 
3, and the positions of the fringes resulting from illuminating the aperture 
of the interferometer with spectrally uniform, nearly collimated light of 
wavelengths, A Q (known wavelength), and A u (unknown wavelength) are shown 
at coordinates x 0 , x u , and xg+i. These coordinates correspond to the fringe 
orders m 0 , m u , and mg+i and wedge heights of h Q , h u , and hg+i- If we 
assume that the medium separating the two surfaces is evacuated (index = 1) 
then the period of the unknown and known fringes is given by 2^ = Au/2B 
and Xq = ^0/23. If the separation of the unknown fringe of order, m u , 
relative to the known fringe of order, m 0 , is given by L = x u - x 0 , then it 
is easy to show that the relationship between these fringes is 

«HjX u = ( m o + L/VAq (^) 

As can be seen from fig. 1, the detected fringes, whose peaks are 
located at the stated coordinates, possess finite widths due to the 
quasi -monochromatic nature of the laser sources and also due to the fringe 
broadening effects of the interferometer (i.e., the instrument finesse). 
Further, since these fringes have been distorted by errors in the detection 
process (e.g., non -uniformity of photosensor response, etc.), the accuracy 
with which the wavelength can be determined depends upon the values of L 
and Xq of the distorted fringes relative to the values which would be 
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obtained if the detection process were perfect. Taking differentials of 
both sides of eq. (1) and assuming that m 0 can be determined exactly by 
calibration procedures with the known laser, one may determine the relative 
uncertainty in the order, m u , to be 

Amy 1 / AL L a2^\ AAq AAy 

= I — + . I + + (2) 

m u m o \ -^o V \> \i 


where we have taken only plus signs as a worst case value. If from other 
means (e.g., a pre-monochrometer, etc.) the order, m u , can be determined 
such that its uncertainty is m u < 0.5, and assuming a^/Xq << AA U /Ay, 
then an upper bound on the relative wavelength uncertainty is given by 


AX, 


1 

m o^o 


AL 



(3) 


From the definitions of L and X 0 the accuracy of this wavelength 
determination depends upon the errors obtained in measuring the various 
coordinate locations of the fringes. 


Estimate of Linecenter by Peak Location 
For a given set of spectral lines, or fringes, which have been detected 
and stored for processing, several methods can be used to measure their 
positions. The particular approach chosen depends upon several factors. 
These include the wavelength accuracy required, any a priori knowledge of 
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the fringe shape, and the time allowed for processing the measurement data. 
One method for this measurement in a system with a diode array readout is by 
peak location. In this approach it is assumed that the peak of the fringe 
irradiance incident upon a perfect array generates the maximum signal 
current. Thus, the determination of the diode location having the maximum 
current locates the peak. As will be shown, however, noise in the detection 
process leads to current fluctuations which distort the true shape of the 
fringe and consequently the peak location. 


Error Analysis 

The error analysis of this method begins with reference to fig. 2 which 
illustrates the fringe irradiance incident on the array. In this figure we 
have neglected the slight asymmetric fringe shape and "ringing" that is 
characteristic of the multiple-beam Fizeau (Ref. 5) and have assumed the 
fringe shape to be Lorentzian, i.e., having the same shape as Fabry-Perot 
fringes (Ref. 6 ). Consequently, we can write the irradiance of a single 
fringe as 


A(x) 



( 4 ) 


where P Q is the optical power contained in the fringe incident on the 
detector plane at x = x Q and 2p is the ful 1 -wi dth-hal f -max (FWHM) of the 
f ri nge. 
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If a sensor element of width. A, is located with its center at x d , and 
has an aperture function given by t(x) = Rect C(x-x d )/A], the optical power 
incident upon the element is the convolution 


P(*d) 


CO . 

C P 0 / 

.1 1 + (x/p 


) 2 


Rect 


x - (x d -x 0 ) 


dx. 


(5) 


A 


Let us assume that x 0 is located within some particular sensor element 
whose center is defined as the center of our coordinate system (i.e., |x 0 | <_ 
0.5 ). Further, let us replace the variable, x, in eq. (5) by x = kA 
where k is continuous and has values at the centers of the sensor elements 
of k d = 0, ±1, +2, ... Using these conditions and defining a diode 
response function, R d (kA), whose value depends upon the coordinate, kA, then 
the current output from the array is 


(k-x 0 /A + 1/2)A 

I(kA) = R d (kA)(P Q /7r) ( dg/p 

J 1 + (S/p) 2 


(k-x 0 /A- 1/2) A 


= R d (kA)(P Q /Tr) 
A/p] I . 


tan" 1 [(k-x Q /A + 1/2) A/p]-tan _1 [( k-x Q /A - 1/2) X 

( 6 ) 


The values at k d = 0, ±1, ±2, give the current from the various detector 
elements. 

If we further define the array response function to be R d (kA) = R 0 + 
r(kA) where R 0 is the response at kA = 0 and r(kA) is the variation of 


M 
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the response about this value, then equation (6) can be expressed as 


I(kA) = (i 0 /e 0 ) 


r(kA) 
1 + R o 


tan 


A/p 

1 + [(k-x 0 /A) 2 -l/4](A/p) 2 


(7) 


In this equation, e 0 = 2tan"^(A/2p) and the current, I 0 , is the current 
obtained when k and x 0 are zero and R d = R 0 (i.e., I 0 = (RoPo0o)/tt ). 

For a perfect array with uniform response, a plot of output current 
from the array for different values of the parameter, |x 0 |/A, is shown in 
fig. 3. Note that for |x 0 |/A = 0.5, the peak current is the same on 
adjacent elements so that the actual location of the peak is uncertain to 
+0.5A. 


The function, R d (kA), present in eq. (7) accounts for those current 
fluctuations which are due to response non-uniformities of the diode array. 
This function can also be used to account for errors which are introduced 
into the measurement by non-uniform spectral irradiance across the Fizeau 
aperture. 

If a small noise modulation is added to the irradiance function given 
by eq. (4) the corresponding optical power incident on a single diode is 
obtained by modifying eq. (5) to include a small power fluctuation, Pn. The 
diode array output current is then 


r(kA) = P(kA) 


R 0 + r(kA) + R 0 . f- 



(8) 


If the last term in this equation is small then the current is 
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I'(kA) - P(kA) 


- 

R 0 + r(kA) + r n ('kA) 


(9) 


where r n = Ro(Pn/P) is the effective response variation in the array due 
to the non-uniform illumination. If we define a total response variation, 
r t = r+r n> ancl ^ et R = ft/ R o» the Previous equation becomes 


I"(kA) * P(kA) R q 1 + R . 


( 10 ) 


In addition to the current fluctuations just described, other 
fluctuations are present in the form of current quantization errors and 
si gnal -to-noise considerations. Specifically, if the maximum signal current 
from the diode array is digitized into Q current intervals of Al = Io/Q, 
then a current error of Aiq = +aI/ 2 is present in the detection process as 
well as a current fluctuation of i a n if si gnal -to-noise effects are 
included. Thus, the signal currents shown in fig. 3 for uniform response 
and no noise must be modified. This results in a peak location uncertainty 
which is greater than our previous result of ±0.5A. Referring to fig. 4, it 
is easy to see that a worst-case combination of these error sources results 
in an uncertainty of the signal peak given by the relation 


I^kA=0; r t (kA) = - r^j -o n - AI/2 = I'^kA; r t (kA) = +rj +a n + AI/2 (11) 


* 
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where I^kA;^) is given by eq. (7) with r replaced by r t . Solving 
this expression for kA gives the coordinate location which satisfies eq. 
(11). The corresponding value of k is 

k = 7 ± 1/2 X 



where we have defined the parameters 1< = x 0 /a, SNR = I 0 /a n > and N = 2p/A 
is the number of diodes that are contained in the fringe width, 2p. Since 
the diode centers are located at integer multiples of A, the value obtained 
in eq. (12) is reduced to the integer which satisfies kj < |k|. 

From the integer obtained above, reference to fig. 4 shows that the 
absolute error, Jax -j J . in locating the peak of a particular fringe at 
coordinate, x-j, is given by Ax -j | = JkjjA - x 0 |. However, since the 
location of x 0 within the k^ = 0 diode is unknown to within ± A/2, the 
value of x 0 must be varied within these limits such that the difference 
| kdA - x 0 is maximum. The maximum absolute error relative to the fringe 
width is thus 



where the error parameter, e, is a function of the wavemeter parameters R, 
SNR, Q, and N. Values of e for various values of these wavemeter 
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parameters are shown in the appendix. The particular value of x which 
maximizes e is also shown. The values for R (0.05, 0.03, and 0.01) were 
chosen as somewhat representati ve of typical devices and employing 
progressive degrees of response calibration (Ref. 7). In addition to the 
table of values in the appendix, a plot of e vs. the number of diodes in 
the width, 2p.j, is shown in fig. 5 for some specific wavemeter parameters. 
Notice that the discrete nature of the diode sampling causes fluctuations in 
the error parameter until the fringe width is spread over about 10 diodes. 

As an example of these results an experimental situation in which SNR - 
100, Q = 256, a total non-uniformity of R = 0.05, and N =10 diodes generates 
a value of e = 0.184. The fringe center can thus be determined to about 
1/5.4 of the FWHM by the fringe peak method. If R is reduced to 0.01 by 
suitable calibration of the diode response and beam non-uniformities are 
negligible this method can determine the fringe center to approximate 
one-tenth of the fringe width. 

Wavelength Error 

The results of the preceding analysis estimate the errors, |Ax^ |, that 
one might expect in locating the various fringe centers by using the fringe 
peak method. In order to translate these results into wavelength errors 
through use of eq. (3), we proceed with some statistical notions regarding 
these coordinate errors. 

As can be seen from fig. 1, the errors, AL and AX^, present in eq. 

(3) are related to the errors in the fringe locations, Ax.,-, by AL = A(x u -x 0 ) 
and AX^ = A(xo+i-x 0 ). If the fringes to be measured are spread across a 
large number of diodes, then repeated measurements of Ax^ would result in 
a histogram or probability density function, p(Ax 1 -), which increases 
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monotonically toward the origin, i.e., it would be peaked at Ax-j = 0. A 
more conservative situation for these measurements would be to consider the 
errors as random. The probability density function would be a rectangular 
distribution where p(Ax-j) = 0 when | Ax -j | > e(2p-j) and P(Ax-j ) = l/e(2P n -) 
when Ax-j | _< e(2p^). This distribution might be reasonable when the number 
of diodes in the FWHM is relatively small and for situations where no 
a priori knowledge of the fringe location is assumed. For such a 
distribution, the mean square value of Ax-j is given by the expected value 
of (Ax.j) 2 , i.e.. 


E[(Ax.j ) 2 ] = J (Ax..) 2 p( Ax . ) d(Ax i ) 


2epi 

*ip t r / 


-2epi 


(Ax.) 2 d(Ax.) = ( 

1 1 \ /T 


/ 


(14) 


The mean-square error obtained in eq. (14) can now be used with eq. (3) 
to find the mean-square error of AA U /X U . Assuming that the coordinate 
error, Ax u , is independent of the errors, Ax Q and AXq+i, but that some 
correlation exists between the various orders of the known laser, an upper 
bound on the rms error of AA U /A U can be shown to be 


AX 


rms 


^ m o\ 


^u 


) 2 + 5 ( 2 p 0 ) 2 


1/2 


(15) 


The widths, 2p^ and 2f^, in this equation are the spatial FWHM 
values. They can be converted to spectral FWHM values, SA U and SA 0 , by 
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the Fizeau relationships which relate the free-spectral -range, order number, 
and spatial period of the fringes at the respective wavelengths, A u and 
A 0 . If the wavelength of the standard or reference laser is chosen 
sufficiently close to the unknown laser wavelength, then eq. (15) becomes 


AX. 


rms 


< 


e ? ? 1/2 

— [5(6x r + (6\ n . 

7K V U 


( 16 ) 


The spectral widths in the previous equation contain both the laser 
linewidths of the reference and unknown laser and also the resolution width 
of the wavemeter. If we assume that reference laser linewidth is small 
compared to the resolution width, and that the fringe shapes are Lorentzian, 
then the absolute rms error is 


AXu 


rms < [5(6X D ) 2 + (SX D + 6A^) 2 ] 


1/2 


A 


( 17 ) 


where <SA r is the Fizeau resolution and 6X^ u is linewidth of the 
unknown laser. Thus, for specific values of A\i|rms» 6A R , and AA £u> the 
required value of e can be determined and the corresponding values of the 
parameters R, SNR, Q, and N found. 

Figure 6 shows an application of eq. (17) in which | AX^ | rms is 
plotted against the error parameter, e, for a laser linewidth of 1.0 pm and 
three values of wavemeter resolution. For the case of an experiment which 
requires absolute knowledge of the wavelength to within 0.25 pm and a 
resolution of 0.50 pm, eq. (17), and fig. 6, indicate that fringe locations 
must be determined to within 0.231 of their respective widths. From the 
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appendix, this value can be obtained by several combinations of the various 
parameters. One such combination which satisfies this value is for a sensor 
where the total non-uniformity of the array and input illumination has value 
the of R = 0.05, and the maximum current from the fringe is digitized by 
an 8-bit analog-to-di gital converter (Q = 256). For these conditions, the 
value of e = 0.231 can be obtained if SNR > 26, and the fringe is spread 
over at least 10 diodes. From fig. 5, any combination of these parameters 
which results in the curve lying below the value of e = 0.231 will achieve 
the wavelength accuracy required in this example. 

Estimate of Linecenter Using Two-Sensor Position Average 

The centers of the fringes which are used to determine the wavelength 
can be found from procedures other than the peak method just described. One 
method which can be used involves searching the array for two diodes (k^* 
and k^ - ) whose output currents are equal to within a specified tolerance. 

The linecenter, x c , is then given by the average position of these two 
diodes, i.e., x c = (k^ + k c j")A/2. 

Analysis of this method for errors in determining the fringe center can 
begin with reference to fig. 3 where the current output from an error-free 
diode array is shown. For the two cases where x 0 = 0 and x 0 = 0.5QA, the 
position average of various diode pairs having equal currents gives a fringe 
center location which corresponds exactly with the true fringe center (x c = 
x 0 ). For x 0 = 0.25A (or any other value between 0 and 0.50 A) no two 
diodes have equal current outputs and clearly some current interval, or 
tolerance must be specified. 

The value of this current interval is somewhat arbitrary but clearly 
the smallest interval that could be chosen is the quantization interval. 
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I 0 /Q , described previously. For actual cases where noise and error sources 
are introduced into the fringes, searching for two diodes whose current 
difference is this small could lead to difficulty. One value of a specified 
current difference which could be used is the relative difference between 
two successive diodes in the FWHM region of the array output. For x 0 - 0 
and for a perfectly uniform array response (R=0), this relative current 
interval can be shown to be given approximately by 1/N where N is the number 
of diodes contained in the FWHM. 

To determine a worst-case error for finding the fringe center we 
proceed in a manner similar to that used in the fringe peak method. We are 
searching for two currents at diode locations x^ 4- = k^A and x^ - = k^~A 
whose outputs due to both the irradiance incident on the array and all the 
various noise sources are equal to within some margin, 5l m . With 
reference to fig. 7, the current output from the array, ^'(k A) (see eq. 
(10)), is shown plotted for an error-free system having a fringe center at 
x 0 = 0. If the various errorswhich were discussed in the previous analysis 
are now introduced into the system, a worst-case error exists by the 
relationship 

V< k d ;\>>- 5I d - si m = + 5I d (18) 


where 
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< - V< k d A * x o> - x 0 . - R ) - a N - AI/2] 

6I d = n d '( k dA; x Q> R) + a N + AI/2] - I d '(k“A; x Q ) 


( 19 ) 


Substituting eq. (7) for I^', equations (18) and (19) combine to give a 
worst-case error condition 



The average location of these two diodes at xt and xT is qiven by 

d d 

x c = ( k t +k T)A/2 resulting in an absolute error of the fringe linecenter of 
d d 

J x o _x c | • absolute error relative to the width, 2p, is defined as in 

eq. (13), then one can show the following relationships: 
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e 



2p 


k d Vx 

P 


o 


z-2e 


k d A ~ X o 

P 


-z-2e 


( 21 ) 


where we have defined the parameter, z, by z = (k^-k^j )/N . With these 
equations, one may solve eq. (21) for the error parameter, e. 


(Hz 2 ) 2 [2 1 6I M 

8z SNR Q I. 



where we have assumed N >> 1 and e/z << 1. 


( 22 ) 


Specific Cases 

The error parameter, e, for this method of determining the fringe 
center can be compared to some values obtained with the fringe-peak method. 
If we assume the current interval, 61^/ I 0 , is given by 1/N as discussed 
previously and the two sampled diodes (k + and k") are at approximately the 
FWHM points, then z 1, and eq. (23) becomes 


e(z=l ) = 1/2(R + 2/SNR + 1/Q + 1/N). 


(23) 
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For Q = 256 and R = 0.05, we have the following table of values of e(z=l) 
vs. SNR: 


Error Parameter vs. SNR and N 


SNR 


e(z= 



N = 10 

N = 20 

N = 30 

N = 40 

5 

0.28 

0.25 

0.24 

0.24 

10 

0.18 

0.15 

0.14 

0.14 

20 

0.13 

0.10 

0.094 

0.089 

30 

0.11 

0.085 

0.077 

0.073 

50 

0.097 

0.072 

0.064 

0.059 

75 

0.090 

0.065 

0.057 

0.053 

100 

0.087 

0.062 

0.054 

0.049 


Comparison of the results in this table against the corresponding 
values shown in the appendix for the fringe-peak method shows a smaller 
error. For example, for N = 10, Q = 256, R = 0.05, and SNR = 100, the 
fringe peak method was shown to give a value of e = 0.184 whereas this 
method gives e = 0.087 or a decrease of about a factor of two in the error 
parameter. However, it must be remembered that we have applied this 
two-position sensor average procedure to relatively smooth fringe shapes. 

In many cases, especially those where the wavelength of pulsed, multi-mode 
lasers is to be determined, the shapes are more complex and often have 
multiple, asymmetric peaks. In addition, the pulse-to-pulse repeatability 
of these complex shapes have not been characterized sufficiently for 
inclusion in this error analysis. For these reasons, then, one would expect 
that the improved accuracy shown in the latter method would decrease and 
approach that obtained by measuring only the peaks. It may be possible, 
however, to regain some of this accuracy by using more than one pair of 
pixel outputs in a statistical averaging process. 
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Conclusions 

The wavelength accuracy that can be obtained from wavemeters employing 
multiple-beam Fizeau interferometers and CCD-type readouts is dependent upon 
several parameters. The non-uniformity of the detector array response and 
fluctuations found in the wavemeter input illumination affect the accuracy 
as well as the current fluctuations in the output due to noise sources such 
as current quantization errors and electronic signal-to-noise ratio. 

Since the wavelength is determined by measuring the locations of fringe sets 
generated by the wavemeter, the method used to determine the fringe 
coordinates also affects the accuracy. The two procedures investigated in 
this paper showed that greater accuracy can be obtained on smooth, 
repeatable laser pulses by the two-position sensor average method, but that 
a simpler technique exists in which only the peak locations are measured. 
Furthermore, when complex, non-repeatable fringe shapes are included, the 
accuracies of the two methods are similar unless a statistical averaging 
process is applied to more than one sensor-pair. 
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Figure 2.- Single fringe irradiance and corresponding 

optical power incident on the diode elements. 








Figure 4.- Error in true peak location due to current fluctuations. 
Actual location of x n is unknown to + 0.50A. 
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Figure 5.- Error parameter vs. number of diodes in fringe width. 
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Figure 7.- Geometry for estimating fringe center using 
two-sensor position average. 
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Error Parameter vs. System Parameters 
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